Current induced Spin Torque in a nanomagnet 
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In a nanomagnet very small polarized currents can lead to magnetic reversal. Treating on the 
same footing the transport and magnetic properties of a nanomagnet connected to magnetic leads 
via tunneling barriers, we derive a closed equation for the time evolution of the magnetization. The 
interplay between Coulomb blockade phenomena and magnetism gives some additional structure to 
the current induced spin torque. In addition to the possibility of stabilizing uniform spin precession 
states, we find that the system is highly hysteretic: up to three different magnetic states can be 
simultaneously stable in one region of the parameter space (magnetic field and bias voltage). 
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Nanomagnets (typically 3 — 4 nm thick grains of mag- 
netic materials with a total spin So ~ 1000) are model 
systems for the study of the interplay between transport 
and magnetism. Their magnetization can be considered 
as a large macrospin without spatial variations 0] (the 
magnetic exchange length is of the order of 7 nm in cobalt 
so that degenerate magnons with non-zero momentum 
can be ignored). Moreover, their mean level spacing is 
of the order of lmeV, bigger than temperatures accessi- 
ble experimentally so that the discreteness of the grain 
spectrum can be resolved Q . A first set of experiments 
used single electron tunneling spectroscopy (connecting 
the grain to electrodes with tunnel junctions) to probe 
the grain's transport properties 0, E| • The differential 
conductance- voltage characteristic of such a device dis- 
plays peaks corresponding to excitations of the grain from 
which the grain's spectrum can be measured 0, 0. In 
a second set of experiments the magnetization of the 
grain is measured directly using a microSQUID technique 
1 . (i . Reversal of magnetic moments as small as 1000 /is 
have been observed providing information on the static 
(switching field as a function of the external field direc- 
tion that allowed the reconstruction of the anisotropy 
tensor) and dynamical (switching times) magnetic prop- 
erties of the system. 

Magnetic reversal by injection of a spin-polarized cur- 
rent was predicted in 1996 by Slonczewski p| and ob- 
served later in magnetic multilayers [8(. A very large in- 
terest for these systems has emerged since, due in partic- 
ular to their large potential for practical applications (in- 
cluding current driven Magnetic Random Access Memo- 
ries and Radio Frequency components) motivating; works 
on magnetic tunnel junctions [qL llOl fllj. The spin dy- 
namics in those systems is usually modelized with the 
Landau-Lifshitz-Gilbert (LLG) equation to which a spin- 
torque term is added 0. 

In this letter, we study the spin dynamics of a nano- 
magnet connected to magnetic leads through tunneling 
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FIG. 1: Nanomagnet connected through tunneling barriers 
to two (magnetic or not) leads with chemical potential and 
magnetization directions fhi (in arbitrary directions). The 
grain magnetization direction is m (||m|| = 1) and its easy 
axis z. The magnetic field H is along z. 



contacts (see Fig. QJ, in the sequential tunneling regime. 
We show that the Coulomb blockade 01 strongly affect 
the spin dynamics in this system. As a consequence, 
the standard LLG equation has to be replaced by Eq. @ 
which is the main result of this letter. In the follow- 
ing, we first outline the derivation of this equation, then 
compare it to the standard LLG equation and study its 
phase diagram as a function of the bias voltage Vt>i as and 
magnetic field H . We find that highly hysteretic current- 
voltage characteristics and metastable magnetic proces- 
sional states could be observed in those systems.. 

Due to the interplay of Coulomb effect and spin dy- 
namics, magnetism and transport have to be treated on 
the same footing in the derivation of the equation govern- 
ing the spin dynamics. We first consider a simple model 
for the isolated nanomagnet, i.e. a mean-field picture of 
the Coulomb blockade, a Stoner instability to get a mag- 
net and a phenomenological modelisation of its intrinsic 
relaxation processes; we then couple the nanomagnet to 
the leads and extract the spin dynamics at lowest order 
in this coupling using a Master equation approach. A 
closely related model has been studied in [13j when the 
exchange interaction is below the Stoner instability. The 
Hamiltonian of the isolated grain has been previously an- 
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alyzed in 0, Q, 0] and reads : 

H g = Y^e a di <T d aa + E(N)-JS-S-^ r S 2 z -h 1 HS z (1) 

bo 

where d^ creates an electron on a one-body level a with 
energy e a , N is the number of electrons in the grain, 
E(N) = Ec(N — N g ) 2 , with Ec the charging energy and 
N g an offset (controlled by a gate voltage), 7 is the gyro- 
magnetic ratio, H the magnetic field. The exchange en- 
ergy (oc J) and uniaxial anisotropy (oc k) depend on the 

total spin S = § Y,a,<T U <n ^ta^i^a^ (<? are Pauli ma- 
trices). Upon increasing J the grain undergoes a Stoner 
instability and acquires a macroscopic spin Sq (equal to 
half of the number of singly occupied one-body states) . 
Its low energy excitations \A) are entirely characterized 
by the set of occupation number n a = 0, 1, 2 of the single 
energy levels together with the spin S z along the easy 
axis jj, fill Hif . Intrinsic relaxation processes (like spin- 
orbit or phonon-magnon scattering) are described phe- 
nomenologically via a weak coupling Tb of the form S + 4> 



to a bosonic bath </> of spectral density pb(E). 

The grain is then coupled via tunneling barriers to 
magnetic leads and we note the tunneling rate of an 
electron with spin a (T or J, along the z axis) coming 
from lead a (Left or Right). For simplicity, we restrict 
ourselves to tunneling events occurring on only one one- 
body state ao, i.e. to one step in the current-voltage 
Coulomb staircase (k B T <C A, where A is the mean level 
spacing) . For larger systems where more levels come into 
play, s pin accumulation effects have to be taken into ac- 
count jlfij . This state can be a majority level (i.e unoc- 
cupied, n ao — 0) or a minority level (n ao = 1). The dif- 
ferent tunneling events (or interactions with the bosonic 
bath) induce a random walk in the Hilbert space of the 
isolated grain, its inner state \A) changing at each event. 
This dynamics is described by the Master equation for 
the probability Pa(^) for the system to be in state \A) at 
time t. The first part is a direct generalization [1J| of the 
standard theory of sequential tunneling |12j . while the 
second part is due to the coupling to the bosonic bath : 
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where AE = Ea — E A > is the energy difference between 
state \A) and \A'), n(E) = 1/(1 + e E ' kT ) (n B (E) = 
\j(e E l kT - 1)) is the Fermi (Bose) function, T the tem- 
perature and p, a the chemical potential of lead a. For 
Eq.(J2J) to be valid, the different tunneling events (and 
interactions with the bosonic bath) have to be incoher- 
ent, which requires T,T B <C kT/h. In addition, Eq.J2J 
does not include the dynamics of the magnetization in 
the xy plane (off-diagonal terms of the density matrix) . 
This is valid when internal precession is fast enough to 
average out contributions transverse to the easy axis. 
The time needed for S z to change significantly is of or- 
der So/T (one tunneling event can only change S z by 
half a unit) so that this condition implies T kSq/%. 
In this paper, we restrict ourselves to that case. The 
effective tunneling rates T„ along the easy axis z can 
be related to the tunneling rates for electrons polar- 
ized along the magnetization of the leads rf,/3 = U,D 
(Up,Down) through = cos 2 ((9 a /2) + ff sin 2 ((9 a /2) 
and T l a = ff cos 2 ((9 Q /2) + sin 2 (6» a /2) where d a is the 
angle between the z axis and the magnetization of lead 
a. The are given by the Fermi golden rule. A com- 
plete derivation of Eq.Q from the microscopic model 



will be provided elsewhere ^t|- Since Sq is macroscopic, 
we proceed with expanding Eq.@ with respect to 1 /So, 
introducing the continuous variable m z = S z /Sq. To ze- 
roth order in 1/So the magnetization does not change and 
the charge equilibrium (between Nq and No + 1 electrons 
in the grain) takes place. Denoting by P+(m z ) (resp. 
P_(m»)) the probability for the grain to contain No + 1 
(resp. No) electrons, it reads: 



n(AE a - p a )P_ 



l-n(AE a -p a ))P 



where AE a (m z ) = e — nam z — fry Ha/ 2 is the energy 
to add an electron in the nanomagnet with spin a. The 
offset energy e includes the charging energy, a possible 
capacitive coupling to an additional gate and the en- 
ergy of the level ao- tj = 1 (resp. -1) when ao is a 
majority (resp. minority) level. On short time scales 
oc h/T, the charge equilibrate and P+/P_ reaches its 
stationary value. The latter is inserted into the first or- 
der part of the 1/Sq expansion to obtain an equation 
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for P(m z ) = P + (m z ) + P_(to z ) of the form dtP(m z ) = 
d mz [-R{m z )P{m z ) + ±- o d mz {D{m z )P(m z )}]. This is a 
Fokker-Planck equation in the Ito form [18|. The asso- 
ciated Langevin equation for m z reads d t m z — R(m z ) + 
^2D(m z )/So £(t), where £(t) is a white noise. The 
noise part is small at large So (the complete expression 
of D(m z ) will be given in |lj|) so that in this letter, 
we focus on the deterministic part R(m z ). Defining the 
asymmetry parameter q = (J2 a 1^ - T a )/(J2 a 1^ + T a ), 
the total tunneling rate T = J2 a rl+T^ and the functions 
fa{mz) = Sa^a n (^-^ <T — we obtain the equation 
governing the spin dynamics of the nanomagnet: 



dm z 
~dt 



(1 - ml) 
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TT 2K 
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where ag oc Tg is the Gilbert damping associated to 
the bosonic bath. Equation (JSJ is the main result of this 
letter. 

Let us now compare Eq. J2J with the LLG equation 
that has been widely used in multilayer systems, 
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It consists of (i) a conservative term (the magnetiza- 
tion precesses around trajectories of constant energy 
E(m) — —nSoml — h^SoHmz); (ii) a phenomenological 
Gilbert damping term a that allows the system to re- 
lax to its equilibrium; (iii) a current induced spin torque 
term Q , written here assuming a magnetic left lead (with 
mi || z) and a non-magnetic right lead. The torque is 
asymmetric and modulated by < q < 1 (q = however 
in a tunneling junction 0, Using the cylindrical 

symmetry of Eq. , we obtain : 



di 
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dt 
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1 + qm z 
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Even though Eq. @ looks similar to Eq.©, its third 
term is more complex leading to richer physics : i) con- 
ducting channels for up/down electrons can open or close 
when m z varies : there is a magnetic blockade induced by 
the anisotropy, and the torque term varies abruptly with 
the bias voltage as the different channels open; ii) this 
term is non-zero even at zero current; it is more than the 
current-induced spin-torque effect, and contains in par- 
ticular relaxations processes due to the leads (13 ■ We 
emphasize that while the usual treatment of spin torque 
physics is done in two separate steps (magnetization dy- 
namics on one hand and transport properties on the other 
hand), here we treat on the same footing the transport 
and magnetic degrees of freedom. Some sterjs in this di- 
rection were taken previously by one of us 14] for grains 
close to equilibrium. 



Let us now investigate the possibility of voltage-driven 
spin reversal. In the framework of the LLG equation, 
a current / with a polarization P flowing through a 
macrospin So reverses the magnetization with a rate 
PI/(eSo) (25o polarized spins +1/2 are needed to fully 
reverse a large spin So from S z = — So to S z = +Sq 
provided the system "absorbs" all the injected polar- 
ized spins). This torque term Pl/eSo competes with 
the relaxation rate an/h (as opposed to k alone in 
magnetic-field-driven reversals) so that the reversal oc- 
curs for 6llg = Plfi/ieSaan) ~ 1. For example, for 
the nanopillars studied in / = Aj (A cross section of 
the pillar and j current density) while So = M s Ad/(h r y) 
(M s saturation magnetization and d thickness of the 
magnetic layer) leads to a critical value of the current 
density j c — eM s dan / (Ph 2 ^) which gives the observed 
value j c ~ 10 7 A cm -2 for typical parameters. In Eq. 
Pjl. the general condition for voltage-driven spin rever- 
sal (as a function Vbias and H) is 9 > 2(1 — q) where 



/(8KCksSo) (correspond- 



e = h[(i- q )r{ + 

ing to the value of the torque when all channels are 
open [l7jL 



The complete phase diagram of Eq. © as a func- 
tion of Vbias and H is found by studying its stable fixed 
points (i.e. the values of m 2 satisfying R(m z ) = with 
d m ^R(m z ) < 0). We restrict the following discussion to 
low temperature, r\ = 1 and to the case fn = (Coulomb 
blockade situation, e is the largest energy). Two identi- 
cal magnetic leads with anti-parallel magnetizations is 
the most favorable case for voltage-driven switching: not 
only up electrons have a higher probability to enter the 
grain than the down electrons, but in addition, once in, 
they have a lower probability to get out so that the in- 
balance between up and down current is enhanced. In 
this case however, the asymmetry parameter q = van- 
ishes. The situation gets more interesting when q ^ 
which is most easily achieved by taking one lead to be 
magnetic while the other is not. The torque can be strong 
enough to destabilize the m z = —1 configuration but 
since its magnitude decreases with increasing m z , a sta- 
ble fixed point —1 < to* < +1 can appear. These fixed 
points are called (uniform) spin precession states (SP) in 
the following. A typical phase diagram is given in Fig. El 
(where the different regions are labeled by the list of the 
stable fixed points). In addition to the regions that can 
exist in the LLG equation, (SP), (SP, +1) and (-1, SP), 
a region (— 1,SP, +1) where three states are stable can 
also exist. In this region, the up spin channel is blocked at 
m z = — 1 so that the torque term does not destabilize the 
(—1) phase. For m z slightly higher, this channels opens 
up, and stabilizes a spin precession fixed point. The con- 
dition for observing a (SP, +1) and (—1, SP, +1) region 
is > (1 — q) 2 /q; the condition for observing a (SP) 
and (-1,SP) region is 6 > (1 - q 2 )/q [13 ■ This phase 
diagram is highly hysteretic as can be inferred from the 
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FIG. 2: (Color online). Phase diagram of Eq. as a func- 
tion of bias voltage Vbias and magnetic field H , measured re- 
spectively in unit of n/e and 2k/"/. The tunneling rates are in 
unit of 8So as and Vbias has been translated of the charg- 
ing energy. ksT = 0.01k. The color code corresponds to the 
stable fixed points as indicated on the figure (e.g. (— 1, +1) 
means that m z — ±1 are stable). We consider one mag- 
netic (left) and one normal (right) lead: Tj cft = 12, = 4, 



^right ~ ri ight = 8. SP denotes a stable Spin Precession state. 
Inset: zoom of the main figure. 



measurement of the current / flowing through the grain. 
In Fig. |3 we show two hysteretic loops corresponding to 
cross-sections of Fig- El We note that in the lower panel, 
we find two different unconnected loops: if the system is 
initially in m z = +1, it will stay there, no matter the 
value of Vbias ■ On the contrary, if the system is initially 
in a SP state or in m z = —1, it will switch, hysteretically 
between the SP and the —1 state. 

In actual nanometric cobalt grains, precise measures 
of the anisotropy K have been obtained using a, vari- 
ety of techniques (including microSQUID Q, \2t\ . tun- 
neling spectroscopy 0, and application of superparam- 
agnetism j2^|) and is of the order of ft/fcs ~ WOmK 
(mainly due to surface anisotropy). Typical currents 
are of a few nA which gives V w 10 10 — 5.10 10 s _1 while 
the total spin So ~ 1000. Typical polarizations are 
P ~ 0.5 — 0.8. More difficult to estimate is the dampin 
constant cub- Switching times were measured in Rcf. 
using a microSQUID technique combined with radio- 
frequency field pulses in the GHz range leading to an 
estimated total damping constant of the order of 0.01 
for a 20nm particle. A significant part of this damp- 
ing is due to spin relaxation through the lead and shall 
not be taken into account in our estimation of the in- 
trinsic damping as (in thin films for instance, the ad- 
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FIG. 3: Current as a function of bias voltage Vbias with 
H = —4.1 (lower panel) and magnetic field H with Vbias = 5.5 
(upper panel) for the same system as Fig. |5] 



junction of a metallic layer below and above a mag- 
netic layer is known to increase the damping by a factor 
10 due to spi n relaxation through the conducting elec- 
trons [20I lijjj). Also, 20 nm particles probably have 
modes of relaxations that disappears for smaller sizes. 
Hence, we estimate as ~ 0.001 — 0.005 (0.005 being the 
bulk value for Cobalt). Putting everything together, we 
get ~ 1. This is of the order of the critical value for 
voltage-driven spin reversal but is probably too small to 
excite spin precession states. Better candidates are mag- 
netic alloys that would allow to decrease the anisotropy 
(permalloy) and/or the total magnetization (weak ferro- 
magnets like NiPd should be good candidates |24|) and 
increase the parameter by one or two orders of mag- 
nitude. Magnetic semiconductors (like GaMnAs) could 
also be material of choice to study the effects described 
in this letter. 

We thank K. Mallick, F. Portier and P. Roche for useful 
comments. 
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